Energy and CO2 fluxes in temperate deciduous broadleaf forests: using local add-covariance data for a global optimization of ORCHIDEE model

Non peer reviewedPublisher PD

Constraining a global ecosystem model with multi-site eddy-covariance data

Peer reviewedPublisher PD

Improving Estimates of Gross Primary Productivity by Assimilating Solar-Induced Fluorescence Satellite Retrievals in a Terrestrial Biosphere Model Using a Process-Based SIF Model

Abstract Over the last few years, solar-induced chlorophyll fluorescence (SIF) observations from space have emerged as a promising resource for evaluating the spatio-temporal distribution of gross primary productivity (GPP) simulated by global terrestrial biosphere models. SIF can be used to improve GPP simulations by optimizing critical model parameters through statistical Bayesian data assimilation techniques. A prerequisite is the availability of a functional link between GPP and SIF in terrestrial biosphere models. Here we present the development of a mechanistic SIF observation operator in the ORCHIDEE (Organizing Carbon and Hydrology In Dynamic Ecosystems) terrestrial biosphere model. It simulates the regulation of photosystem II fluorescence quantum yield at the leaf level thanks to a novel parameterization of non-photochemical quenching as a function of temperature, photosynthetically active radiation, and normalized quantum yield of photochemistry. It emulates the radiative transfer of chlorophyll fluorescence to the top of the canopy using a parametric simplification of the SCOPE (Soil Canopy Observation Photosynthesis Energy) model. We assimilate two years of monthly OCO-2 (Orbiting Carbon Observatory-2) SIF product at 0.5° (2015?2016) to optimize ORCHIDEE photosynthesis and phenological parameters over an ensemble of grid points for all plant functional types. The impact on the simulated GPP is considerable with a large decrease of the global scale budget by 28 GtC/year over the period 1990?2009. The optimized GPP budget (134/136 GtC/year over 1990?2009/2001?2009) remarkably agrees with independent GPP estimates, FLUXSAT (137 GtC/year over 2001?2009) in particular and FLUXCOM (121 GtC/year over 1990?2009). Our results also suggest a biome dependency of the SIF-GPP relationship that needs to be improved for some plant functional types.Peer reviewe

Land surface model parameter optimisation using in situ flux data: comparison of gradient-based versus random search algorithms (a case study using ORCHIDEE v1.9.5.2)

Land surface models (LSMs), which form the land component of earth system models, rely on numerous processes for describing carbon, water and energy budgets, often associated with highly uncertain parameters. Data assimilation (DA) is a useful approach for optimising the most critical parameters in order to improve model accuracy and refine future climate predictions. In this study, we compare two different DA methods for optimising the parameters of seven plant functional types (PFTs) of the ORCHIDEE LSM using daily averaged eddy-covariance observations of net ecosystem exchange and latent heat flux at 78 sites across the globe. We perform a technical investigation of two classes of minimisation methods – local gradient-based (the L-BFGS-B algorithm, limited memory Broyden–Fletcher–Goldfarb–Shanno algorithm with bound constraints) and global random search (the genetic algorithm) – by evaluating their relative performance in terms of the model–data fit and the difference in retrieved parameter values. We examine the performance of each method for two cases: when optimising parameters at each site independently (“single-site” approach) and when simultaneously optimising the model at all sites for a given PFT using a common set of parameters (“multi-site” approach). We find that for the single site case the random search algorithm results in lower values of the cost function (i.e. lower model–data root mean square differences) than the gradient-based method; the difference between the two methods is smaller for the multi-site optimisation due to a smoothing of the cost function shape with a greater number of observations. The spread of the cost function, when performing the same tests with 16 random first-guess parameters, is much larger with the gradient-based method, due to the higher likelihood of being trapped in local minima. When using pseudo-observation tests, the genetic algorithm results in a closer approximation of the true posterior parameter value in the L-BFGS-B algorithm. We demonstrate the advantages and challenges of different DA techniques and provide some advice on using it for the LSM parameter optimisation.</p

A review of applying second-generation wavelets for noise removal from remote sensing data.

The processing of remotely sensed data includes compression, noise reduction, classification, feature extraction, change detection and any improvement associated with the problems at hand. In the literature, wavelet methods have been widely used for analysing remote sensing images and signals. The second-generation of wavelets, which is designed based on a method called the lifting scheme, is almost a new version of wavelets, and its application in the remote sensing field is fresh. Although first-generation wavelets have been proven to offer effective techniques for processing remotely sensed data, second-generation wavelets are more efficient in some respects, as will be discussed later. The aim of this review paper is to examine all existing studies in the literature related to applying second-generation wavelets for denoising remote sensing data. However, to make a better understanding of the application of wavelet-based denoising methods for remote sensing data, some studies that apply first-generation wavelets are also presented. In the part of hyperspectral data, there is a focus on noise removal from vegetation spectrum

Quantifying Vegetation Biophysical Variables from Imaging Spectroscopy Data: A Review on Retrieval Methods

An unprecedented spectroscopic data stream will soon become available with forthcoming Earth-observing satellite missions equipped with imaging spectroradiometers. This data stream will open up a vast array of opportunities to quantify a diversity of biochemical and structural vegetation properties. The processing requirements for such large data streams require reliable retrieval techniques enabling the spatiotemporally explicit quantification of biophysical variables. With the aim of preparing for this new era of Earth observation, this review summarizes the state-of-the-art retrieval methods that have been applied in experimental imaging spectroscopy studies inferring all kinds of vegetation biophysical variables. Identified retrieval methods are categorized into: (1) parametric regression, including vegetation indices, shape indices and spectral transformations; (2) nonparametric regression, including linear and nonlinear machine learning regression algorithms; (3) physically based, including inversion of radiative transfer models (RTMs) using numerical optimization and look-up table approaches; and (4) hybrid regression methods, which combine RTM simulations with machine learning regression methods. For each of these categories, an overview of widely applied methods with application to mapping vegetation properties is given. In view of processing imaging spectroscopy data, a critical aspect involves the challenge of dealing with spectral multicollinearity. The ability to provide robust estimates, retrieval uncertainties and acceptable retrieval processing speed are other important aspects in view of operational processing. Recommendations towards new-generation spectroscopy-based processing chains for operational production of biophysical variables are given

Interannual vegetation phenology estimates from global AVHRR measurements - Comparison with in situ data and applications

International audienceIn a previous paper, we described a procedure to correct the directional effects in AVHRR reflectance time series. The corrected measurements show much less high frequency variability than their original counterparts, which makes them suitable to study vegetation dynamics. The time series are used here to estimate the start and ending dates of the growing season for 18 years from 1982 to 1999. We focus on the interannual variations of these phenological parameters.A database of in situ phenology observations is used to quantify the accuracy of the satellite-based estimates. Although based on a limited sampling of the Northern mid and high latitudes, the comparison indicates that i) the satellite phenological product contains meaningful information on interannual onset anomalies; ii) there is a higher degree of consistency over regions covered by Broadleaf Forests, Grasses and cereal Crops than over those covered by Needleleaf Forests or Savannas; and iii) the satellite phenological product is of lower quality in regions with mountainous terrain. In favorable conditions, interannual variations of the onset are captured with an accuracy of a few days.As this satellite-derived product captures the interannual onset variability at ground-truth sites, we confidently use it to larger scales studies. Mapped at a continental scale, the onset anomalies show coherent patterns at the regional (≈ 1000 km) scale for the mid and high latitudes of the Northern hemisphere, which is consistent with a meteorological forcing. In the tropics, there is larger spatial heterogeneity, which suggests more complex controls of the phenology. The relation between vegetation phenology and climate is further investigated over Europe by comparing the variability of the satellite-derived vegetation onset and that of the winter North Atlantic Oscillation index, at a fine spatial scale. The strong correlations observed confirm that this climate forcing parameter explains most of the onset variability over a large fraction of Northern Europe (earlier onsets for positive winter NAO), with lower impact towards the south and opposite effects around the Mediterranean basin. The NAO has a predictive character as the January–February NAO index is strongly correlated with the vegetation onset that occurs around April in Northern Europe

Ecosystem model optimization using in situ flux observations: Benefit of Monte Carlo versus variational schemes and analyses of the year-to-year model performances

Terrestrial ecosystem models can provide major insights into the responses of Earth's ecosystems to environmental changes and rising levels of atmospheric CO2. To achieve this goal, biosphere models need mechanistic formulations of the processes that drive the ecosystem functioning from diurnal to decadal timescales. However, the subsequent complexity of model equations is associated with unknown or poorly calibrated parameters that limit the accuracy of long-term simulations of carbon or water fluxes and their interannual variations. In this study, we develop a data assimilation framework to constrain the parameters of a mechanistic land surface model (ORCHIDEE) with eddy-covariance observations of CO2 and latent heat fluxes made during the years 2001-2004 at the temperate beech forest site of Hesse, in eastern France. As a first technical issue, we show that for a complex process-based model such as ORCHIDEE with many (28) parameters to be retrieved, a Monte Carlo approach (genetic algorithm, GA) provides more reliable optimal parameter values than a gradient-based minimization algorithm (variational scheme). The GA allows the global minimum to be found more efficiently, whilst the variational scheme often provides values relative to local minima. The ORCHIDEE model is then optimized for each year, and for the whole 2001-2004 period. We first find that a reduced (< 10) set of parameters can be tightly constrained by the eddy-covariance observations, with a typical error reduction of 90 %.We then show that including contrasted weather regimes (dry in 2003 and wet in 2002) is necessary to optimize a few specific parameters (like the temperature dependence of the photosynthetic activity). Furthermore, we find that parameters inverted from 4 years of flux measurements are successful at enhancing the model fit to the data on several timescales (from monthly to interannual), resulting in a typical modeling efficiency of 92% over the 2001-2004 period (Nash-Sutcliffe coefficient). This suggests that ORCHIDEE is able robustly to predict, after optimization, the fluxes of CO2 and the latent heat of a specific temperate beech forest (Hesse site). Finally, it is shown that using only 1 year of data does not produce robust enough optimized parameter sets in order to simulate properly the year-to-year flux variability. This emphasizes the need to assimilate data over several years, including contrasted weather regimes, to improve the simulated flux interannual variability. © Author(s) 2014